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ABSTRACT 

We report on a statistical study of the 51 radio galaxies at the millijansky flux level from the Faint 
Images of the Radio Sky at Twenty centimeters, including their optical morphologies and structure 
obtained with the Hubble Space Telescope. Our optical imaging is significantly deeper (~2 mag) than 
previous studies with the superior angular resolution of space-based imaging. We that find 8/51 (16%) 
of the radio sources have no optically identifiable counterpart to AB ~ 24 mag. For the remaining 43 
sources, only 25 are sufficiently resolved in the HST images to reliably assign a visual classification: 
15 (60%) are elliptical galaxies, 2 (8%) are late-type spiral galaxies, 1 (4%) is an SO, 3 (12%) are 
point-like objects (quasars), and 4 (16%) are merger systems. We find a similar distribution of optical 
types with measurements of the Sersic index. The optical magnitude distribution of these galaxies 
peaks at 7^20.7 ± 0.5 AB mag, which is ~3 mag brighter than the depth of our typical HST field 
and is thus not due to the WFPC2 detection limit. This supports the luminosity-dependent density 
evolutionary model, where the majority of faint radio galaxies typically have L*-optical luminosities 
and a median redshift of z ~ 0.8 with a relatively abrupt redshift cut-off at z > 2. We discuss our 
results in the context of the evolution of elliptical galaxies and active galactic nuclei. 
Subject headings: galaxies: starburst — radio continuum: galaxies 



1. INTRODUCTION 

Faint radio galaxies are good probes of typical galaxies 
at large distances since they lack a significant contribu- 
tion from non-thermal light or strong emission lines in 
their optical spectra (eg. Kron, Koo, & Windhorst 1985; 
Benn et al. 1993; Hammer et al. 1995a). Furthermore, 
their selection is relatively unaffected by dust absorption 
when these galaxies are identified by their synchrotron 
emission. In the local Universe, powerful radio sources 
are primarily found in giant elliptical galaxies, which sug- 
gests that radio selection may be an efficient method of 
identifying high-redshift ellipticals. However, at lower 
radio powers, typically below the break in the radio lu- 
minosity function, there is an increasing fraction of star- 
forming and active galaxies (such as Seyferts and LIN- 
ERS; eg. Condon 1989; Benn et al. 1993; Lowenthal 1997; 
Hopkins et al. 2003; Best et al. 2005; Afonso et al. 2006; 
Mainieri et al. 2008). Often of the galaxy identification is 
not based on deep, high-resolution optical imaging, but 
rather on broadband optical colors. Such classification 
schemes are only reliable to the extent to which these 
colors directly map onto the Hubble sequence at a given 
redshift. Therefore, in this paper we focus on the optical 
morphologies and structure of radio galaxies with fluxes 
of 1 mJy< Si a GHz ;$ 100 mJy (we will denote Si a as 
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integrated radio flux at 1.4 GHz throughout) obtained 
from the Hubble Space Telescope (HST). 

Early studies indicate that the majority of millijan- 
sky radio galaxies typically have extended radio emis- 
sion with optical colors and luminosities comparable to 
local elliptical galaxies (Kron, Koo, & Windhorst 1985; 
Windhorst et al. 1985; Benn et al. 1993; Hammer et al. 
1995a). The remainder of these radio sources are gen- 
erally quasars or very blue galaxies, which are often in- 
terpreted as late-type, starburst systems. Like the most 
powerful radio sources, millijansky galaxies have a broad 
redshift distribution of 0<z<1.5 (Overzier et al. 2003), 
which reflect the increase of active galaxies and radio 
emission from z ~ to z ~ 2 (eg. Hopkins & Beacom 
2006). The space density of millijansky radio galax- 
ies evolves strongly with redshift, which can be repro- 
duced by a simple evolutionary model. Condon (1989) 
argue that the majority of these radio galaxies reside in 
a "shell" of radius z ~ 0.8 with thickness of Sz ~ 0.8. 
In addition to supporting this model, Waddington et al. 
(2001) find a deficit of high redshift (z > 2) millijansky 
radio galaxies in the Hercules field of the Leiden-Berkley 
Deep Survey (LBDS). They interpret this result as an 
effective redshift "cut-off," which varies from z ~ 2 for 
the most luminous radio sources, to z ~ 1 for the lower 
luminosity galaxies. 

Millijansky radio sources are intrinsically rather lu- 
minous at optical wavelengths, with absolute magni- 
tudes in a relatively narrow range (±0.8 mag) centered 
at M r * ~ —22 mag and M r * ~ —23 mag for blue and 
red galaxies, respectively (Machalski & Godlowski 2000; 
Ivezic et al. 2002). From photographic plate and deep 
CCD imaging, Windhorst et al. (1984b) and Wadding- 
ton et al. (2000) found the that the apparent magnitude 
distribution from the LBDS peaks around I ~ 22 mag, 
which at z = 0.8 corresponds to Mb — —21.5 mag. While 
there is a good deal of deep ground-based imaging on mil- 
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lijansky radio sources, there is little high-resolution HST 
imaging. Early works with HST are generally based on 
at most a few objects (Pascarelle et al. 1996; Waddington 
et al. 1999, 2002). 

In the microjansky flux range, Afonso et al. (2006) 
identified 64 faint radio galaxies of which 57 had optical 
counterparts in the Advanced Camera for Surveys (ACS) 
imaging of the Great Observatories Origins Deep Survey, 
South (GOODS-S; Giavalisco et al. 2004). The median 
photometric redshift of the optically identified microjan- 
sky radio sources is also z^O.S. Since the redshift distri- 
bution of radio galaxies changes only very slowly for ra- 
dio fluxes below Si.4~l Jy (Windhorst et al. 1990), this 
median photometric redshift is consistent with the evo- 
lutionary model of Condon (1989). Based on their X-ray 
properties and limited optical spectroscopy, the major- 
ity of these microjansky radio galaxies are either active 
galactic nuclei (AGN) or star-forming galaxies (Afonso et 
al. 2006). Furthermore, the optical counterparts gener- 
ally have disturbed morphologies, which suggests many 
are in the later stages of merging. 

In this paper, we investigate the optical properties of 
millijansky radio galaxies selected from the Faint Images 
of the Radio Sky at Twenty centimeters (FIRST; Becker 
et al. 1995). The FIRST survey covers ~ 10 4 square de- 
grees of the north Galactic cap to a 95% completeness 
limit of ~ 2.8 mJy. In order to reliably determine the 
morphology of these small and faint optical counterparts, 
we require the angular resolution of space-based imag- 
ing. Therefore, our sample is derived from the overlap of 
FIRST sources with the existing HST-WFPC2 Archive 
as of 2002. All magnitudes quoted herein are in the AB 
system (Oke & Gunn 1983). This paper is organized 
as follows: In § 2 we outline the critical steps to define 
our final millijansky sample, in § 3 we discuss our four 
primary results, in § 4 we provide general comments on 
individual objects, and in § 5 we conclude with a general 
summary and additional discussion of our results. 

2. PROCEDURE 

2.1. Data Selection 

We begin by correlating the positions from the FIRST 
catalog with the HST-WFPC2 Archive. Out of > 
10 4 WFPC2 images, we obtain ~ 850 potential frames 
with FIRST sources. These WFPC2 images were stacked 
and cosmic-ray cleaned following the method discussed 
by Cohen et al. (2003) to yield a sample of 177 opti- 
cal candidates. The HST-WFPC2 dataset contains all 
the archival F606W (V) and F814W (/) images with ex- 
posure times longer than 160 seconds, however only a 
single band is generally available for a given target, as 
the primary goals of the HST images. Roughly 50% of 
the WFPC2 data are observed as HST parallel fields, 
which are images taken by WFPC2 when a different in- 
strument was observing the primary target. In general, 
these millijansky radio sources are extremely difficult to 
image since they are typically optically faint (/ > 20 mag) 
and distributed over degree scales. Therefore, these par- 
allel observations provide an efficient method to optically 
image this elusive population. 

2.2. Secondary Selection Criteria 

Since the motivation of this work is to study the op- 
tical morphologies and structure of a well defined, ran- 



TABLE 1 

Sample Reduction Summary 



Reason 



Number 



Radio Sources In Final Sample 

Excluded because of Pointed HST Observations 

Excluded because NGC or Nearby Contamination 

Excluded because Incomplete in Radio Flux Cutoff 

Misc a 



51 
43 
32 
47 
4 



Total 



177 



a This category includes fields where the target is uncertain due 
to a galaxy group or cluster, saturated WFPC2 fields, and/or local 
WFPC2 defects precluding inclusion in the optical sample. 

dom, and complete sample of weak radio galaxies with 
the superior resolution of HST, many of the candidates 
must be eliminated from the sample. There are four 
reasons what we could exclude a given field: (1) Major 
WFPC2 image defects, such the majority of the pixels 
are saturated from nearby stars; (2) a significant frac- 
tion of the field-of-view is covered by a nearby object, 
such as an NGC-type galaxy, which would prevent reli- 
able optical flux measurements and morphological classi- 
fications; (3) the WFPC2 image was centered on a dense 
star-field or galaxy cluster to eliminate ambiguity in the 
optical identifications; or (4) the radio source was tar- 
geted by HST. In this way, we can subselect a complete 
sample of background objects from a non-random set of 
foreground observations. Furthermore, owing to the rel- 
atively small WFPC2 field-of-view, most FIRST radio 
sources are generally the primary target. In Table 1, we 
summarize the eliminated fields. The selection effects 
of the various HST programs is not always possible to 
quantify based on the HST target lists and PI proposals 
available in the STScI/HST Archive. In Table 2, we list 
the targeted radio sources that were removed to reduce 
the contamination of preselected objects in constructing 
our random sample. These criteria provide a compro- 
mise between the resulting completeness and reliability 
of the final sample. The majority of the radio sources are 
brighter than the 95% completeness limit of the FIRST 
survey of S 1A ^ 2.8 mJy. From the 177 FIRST candi- 
dates, our final sample contains 51 radio galaxies with 
HST-WFPC2 observations. 

2.3. Astrometric Correction of HST Fields 

Due to technical issues in the the world coordinate sys- 
tem (WCS) generated in the HST-WFPC2 pipeline, the 
astrometry keywords in the FITS headers of WFPC2 
fields prior to 1997 September 15 may be unreliable 
(Biretta et al. 2000). To better calibrate the astrometry 
of these frames, we match the objects in each field with 
the United States Naval Observatory (USNO) A2.0 cata- 
log (Monet et al. 1996), which is the standard in the FK5 
system. We derive shifts in Right Ascension (Aa) and 
Declination (AS) from the positional differences between 
the HST-WFPC2 images and USNO catalog, which we 
show in Figure 1. These offsets are individually applied 
to their respective WFPC2 frame. However, several fields 
lacked usable USNO sources, and were left uncorrected. 
In general, these were high-latitude fields or contained 
only saturated stars, which made accurate determina- 
tion of their centers unreliable. While the dispersion in 
AS is not significantly higher for the pre- 1997 data, it 
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TABLE 2 

Targeted Sources Removed from Sample 



Radio Source a 


RA J2000 


Dec J2000 


log #1.4 


WFPC2 Target 




(deg) 


(deg) 


(mjy) 


(targname) b 


J001303.7-012305 


3.265625 


-1.384777 


0.61 


CL-ULIR00105-0139 


J022448. 1-000711 


36.200542 


-0.120111 


0.04 


GAL-CLUS-022448-000717 


J074123.0+320810 


115.345963 


32.136360 


0.60 


QDOTULIR073810+321511 


J082209.5+470552 


125.539925 


47.098114 


1.72 


A646 


J082313.7+275120 


125.807213 


27.856083 


1.03 


QDOTULIR082010+280119 


J083803.6+505509 


129.515091 


50.919613 


0.88 


QDOTULIR0838+5055 


J090850.9+374818 


137.211304 


37.805500 


3.06 


3C217 


J094308.4+465700 


145.785828 


46.950359 


0.20 


GAL-CLUS-093942+47 


J095254.3+435348 


148.226196 


43.896782 


0.85 


GAL-CLUS-0952+44 


J095644. 7+493034 


149.187088 


49.513168 


0.68 


A895 


J102326.2+124855 


155.859283 


12.815556 


0.22 


ABELL999-BCG 


J103141. 7+350231 


157.924042 


35.044613 


1.64 


A1033 


J105648.8-033727 


164.203873 


-3.623805 


0.82 


GAL-CLUS-105427-032125-POS1 


J105659.5-033727 


164.248123 


-3.624416 


1.09 


GAL-CLUS-105427-032125-POS5 


J110050.3+104654 


165.211929 


10.783639 


1.06 


AO1058+11 


J115106.8+550443 


177.778549 


55.078945 


0.95 


UGC06823 


J122906.7+020308 


187.277954 


2.052416 


4.57 


PG1226+023 


J123357.1+074206 


188.488297 


7.701722 


0.87 


NGC4526 


J125929. 2+361713 


194.875290 


36.284283 


2.17 


6C1257+36 


J131128.9-012116 


197.870880 


-1.354444 


0.36 


ABELL1689 


J131130. 0-012028 


197.875122 


-1.341333 


0.94 


A1689-10 


J131131. 5-011931 


197.880920 


-1.325666 


1.63 


ABELL1689 


J131132. 6-011959 


197.886169 


-1.333000 


0.39 


A1689-10 


J133238.3+503335 


203.159988 


50.560001 


0.78 


A1758G7 


J133239.4+503432 


203.164673 


50.575615 


0.77 


A1758G7 


J133557.2+544338 


203.987244 


54.726387 


1.09 


MS1333.9+5500 


J134337. 1-001525 


205.904648 


-0.257166 


0.74 


GAL-CLUS-134339-001349 


J134447.4+555410 


206.197723 


55.903053 


0.10 


IR13428+5608.0 


J134733. 3+121724 


206.889297 


12.290195 


3.69 


PKS1345+12-PSF 


J135609.9+290535 


209.041580 


29.093473 


1.03 


NGWULIR1353+2920 


J140318.1+542157 


210.824783 


54.365974 


1.14 


NGC5457-FLD1 


J141721. 3+132429 


214.339127 


13.408417 


0.37 


MS1414.9+1337 


J142357.6+383247 


215.990982 


38.546196 


1.02 


FSC142118+3845 


J142553.5+374805 


216.473129 


37.801472 


0.27 


A1914 


J143242.8+245614 


218.172867 


24.921638 


2.09 


B2+1430+25 


J144516.4+095836 


221.318619 


9.976666 


3.38 


1442+101 


J151002.9+570243 


227.512466 


57.045555 


2.41 


GB1508+5714 


J160218. 2+155912 


240.576035 


15.986777 


0.87 


ABELL2147-BCG 


J162124.6+381008 


245.353226 


38.169422 


0.75 


R1621+38 


J162439. 6+234524 


246.162979 


23.753334 


3.15 


GAL-CLUS-3C336-POS2 


J162548.7+264658 


246.453186 


26.782833 


0.99 


Q1623+268 


J164658. 9+454824 


251.745468 


45.806778 


0.74 


QDOTULIR1645+4553 


J231713.1-110034 


349.305084 


-11.009694 


0.82 


CL-ULIR23146-1117 



a The name of the source as per IAU recommendations, as used 
by FIRST. FIRST Jhhmmss.s+ddmmss in which the coordinates 
are equinox J2000.0 and are truncated (Becker et al. 1995). 

b This is the header keyword for the field target from the WFPC2 
images. 

does increase by a factor of ~4 for the dispersion in Aa. 
The applied shifts roughly follow a normal distribution 
centered on ((Aa,A5)) = (-O'.'l, O'/O) with standard de- 
viations of (< pt ,a° pt ) = (0"85,0'/55). This mean offset 



gives a systematic uncertainty of 
random uncertainty is ~l'/0. 



2.4. 



O'/l, while the la 



WFPC2 Image Preparation and Investigation 

We excise 300x300 pixel "postage stamp" images from 
the astrometrically corrected and stacked WFPC2 mo- 
saics centered on each FIRST source position. In Fig- 
ure 2, we show the 51 WFPC2 stamps with FIRST radio 
contours placed at 2 n multiples of 0.31 mjy, and list 
the relevant photometric and morphological information 



in the upper-left corner. Optical counterparts to the 
given radio sources are identified by the spatial align- 
ment between the optical image and the center of the 
radio contours. Since several objects have complex radio 
morphologies (such as elongated or double isophotes), 
the identification of each most likely optical counterpart 
was for such objects done on a case- by-case basis. For 
the FIRST sources with a double-lobed morphology, we 
adopt the optical object which is nearest to the geomet- 
ric center of the lobes as the most likely identification 
(eg. Mackay 1971). If the optical counterpart is consid- 
erably off the line which connects the two lobes, then we 
use the the orientation and outer isophotes to aide in the 
identification (such as in the case of J094930. 7+295938). 
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-1 o 1 
Ac* (arcsec) 

Fig. 1. — Optical astrometric corrections in RA and Dec. 
These relative offsets are computed as the positional difference be- 
tween the HST-WFPC2 and USNO coordinates for multiple stars 
in a given field. Both distributions are roughly Gaussian with 
mean values of (Aa, AS) = (— O'/l, 0"0), with standard deviations 
(cr° pt ,cr° pt ) = (0'/85, 0'/55). The ellipses represent the la, 2a, and 
3a contours. These systematic offsets and uncertainties are used to 
determine the most likely optical identification for radio galaxies 
discussed in § 2.4. 



TABLE 3 

Unidentified Radio Sources 



Radio Source a 


Passband 


Limit b 


log5i.4 






(mag) 


(mjy) 


J030237.6+000818 


F814W 


24.5 


0.70 


J082828.4+344131 


F814W 


24.7 


1.01 


J120326.4+443635 


F606W 


25.0 


1.43 


J121839.7+295325 


F606W 


24.8 


1.53 


J125635.3+215632 


F814W 


24.2 


0.52 


J150109. 7+225111 


F814W 


23.6 


0.82 


J151057.4+312722 


F606W 


23.7 


0.49 


J164911. 4+305226 


F606W 


24.3 


0.44 



a The name of the source as per IAU recommendations, as used 
by FIRST. FIRST Jhhmmss.s+ddmmss in which the coordinates 
are equinox J2000.0 and are truncated (Becker et al. 1995). 

b The limiting AB magnitude in specified filter given as 50% com- 
pleteness limit following Snigula et al. (2002). 

Both lobes of double radio sources are confirmed to not 
separately contain optical counterparts. Stamps with 
no optical object within the radio contours (complete to 
2.8 mjy) were deemed unidentified radio sources (here- 
after, Unid) to the limiting optical magnitude cataloged 
in Table 3. For identified optical counterparts, we assign 
visual types and measure their surface brightness profiles 
(see § 3.1), however the faintest and smallest HST ob- 
jects are generally not classifiable (Odewahn et al. 1996; 
Cohen et al. 2003). 

To quantify the reliability of the radio-optical identi- 
fication, we use the likelihood ratio (LR) technique de- 
scribed by de Ruiter et al. (1977). If we assume that the 
optical and radio sources are intrinsically located at the 
same position, and that potential clustering does not lead 
to a large fraction of false identifications, then the likeli- 
hood that a radio source and its potential identification 



are the same physical object can be defined as: 
where 

/ / Q/rad Q/°Pt \ / (5 ra d (5°Pt x 



0~8 



\ = 7ra a (T5PBG, 



^ = \/K ad ) 2 + K pt ) 2 - 



(1) 

(2) 
(3) 
(4) 
(5) 



where pbg is the surface density of objects, (cr^ ad ,cr£ ad ) 
are the radio positional uncertainties (Becker et al. 1995), 
(o-^ pt ,cr^ pt ) are tne uncertainties in the WCS shifts de- 
rived from the fits headers (discussed in § 2.3 and Fig- 
ure 1), and (a rad ,(5 rad ) and (a opt ,£ opt ) are the celestial 
coordinates of the radio and optical sources, respectively. 

In these relatively shallow and narrow HST fields, an 
accurate background density of objects cannot always be 
determined. Therefore, we adopt the method of Snigula 
et al. (2002) to estimate the galaxy flux at the 50% com- 
pleteness limit. This flux is found by assuming a given 
light-profile, and using the noise properties of the image 
to determine the brightness limit to which 50% the ob- 
jects would have been recovered. For a given total mag- 
nitude, pure disk galaxies have a lower surface brightness 
than pure ellipticals, and are therefore somewhat more 
difficult to detect. However, the 50% completeness limit 
for an exponential and de Vaucouleurs typically differ by 
only < 0.1 mag for a fixed size of r = 0'/3 (Cohen et al. 
2003). Given this minor difference in brightness for the 
two light-profiles, we adopt the limits derived from the 
de Vaucouleurs profile, which are given in Table 4. Fi- 
nally, we integrate the observed /-band number counts 
(Metcalfe et al. 1996; Odewahn et al. 1996; Casertano et 
al. 2000; Gardner & Satyapal 2000; Yasuda et al. 2001) to 
these derived limits to determine the background density 
(Pbg) in equation (3). 

The resulting likelihood ratios are listed in Table 4, 
and their distribution is plotted in Figure 3. Following 
de Ruiter et al. (1977), a LR>2 provides a compromise 
between reliability and completeness for the radio-optical 
identifications. Therefore, we adopt the de Ruiter et al. 
(1977) likelihood threshold, where objects with LR > 2 
are a positive identification. For radio sources with no 
obvious counterpart (the Unids), the likelihood ratio was 
computed for the nearest optical source, and all of these 
systems have a Li?<C2. 



3. RESULTS 



3.1. 



Optical Morphologies 

Many of the counterparts to these faint radio galax- 
ies are also very faint in the optical. Only 2/43 galaxies 
have / < 16.0 mag, the approximate brightness limit to 
classify galaxies from ground-based imaging in average 
seeing conditions. Of the 43 radio sources for which an 
optical counterpart was found, only 25 can be visually 
classified to a brightness limit of AB < 22 mag. The vi- 
sual morphologies were determined by three independent 
observers (J.R., S.H.C., and R.A.W.), and the average 
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TABLE 4 


















Flux Results 










Radio Source a 


Band b 


Magnitude 


SB C 


logSi. 4 d 


Depth 6 


LR f 


i's 


W ~ r'Y 


(r' - i'Y 






(mag) 


(mag/D") 


(mjy) 


(mag) 




(mag) 


(mag) 


(mag) 


J002219.1-013030 


V 


24.59 


23.39 


0.47 


25.2 


11.17 








J004322. 3-001343 


V 


18.93 


16.12 


0.83 


23.9 


10.17 


18.91+0.02 


0.28+0.02 


0.12+0.03 


J012616. 6-012126 


V 


19.84 


22.45 


1.15 


23.9 


2.03 








J030237.6+000818 


I 






0.70 


24.5 


0.0 








J030249.5+000615 


I 


19.79 


20.33 


0.15 


24.7 


6.42 


20.13+0.18 


0.64+0.27 


0.65+0.22 


J082820. 6+344321 


I 


25.59 


27.21 


0.69 


24.9 


11.65 








J082828.4+344131 


I 






1.01 


24.7 


0.0 








J084715.5+443752 


I 


24.66 


23.27 


0.62 


23.5 


44.39 








J084849.5+445550 


I 


19.29 


21.06 


0.67 


24.6 


10.53 


19.13+0.06 


2.40+0.61 


1.02+0.12 


J091205.2+350506 


V 


23.79 


22.50 


0.93 


23.7 


17.04 








J091251.0+525928 


I 


21.14 


20.80 


0.47 


22.9 


23.59 


20.42+0.12 


4.79+0.89 


1.51+0.34 


J094926.6+295941 


V 


20.06 


21.85 


1.25 


25.4 


8.43 


19.14+0.07 


2.08+0.28 


0.64+0.12 


J094930.7+295938 


V 


22.70 


22.54 


0.55 


25.6 


3.87 








J100354.5+285911 


V 


23.10 


24.77 


0.60 


24.6 


6.25 


22.44+0.44 


0.96+1.45 


1.70+1.07 


J102437.2+470312 


I 


21.75 


20.36 


0.39 


24.3 


15.62 


22.31+0.31 


1.36+0.78 


0.91+0.49 


J102744.6+282921 


I 


22.44 


24.31 


0.15 


24.6 


13.33 








J103452.3+394704 


I 


21.73 


21.94 


1.40 


24.5 


27.44 


21.82+0.48 


1.92+1.26 


0.39+0.63 


J104630.8-001213 


I 


18.8 


21.73 


0.65 


24.0 


28.24 


19.19+0.08 


1.91+0.29 


0.71+0.13 


J104757.0+123835 


V 


20.61 


20.77 


0.94 


24.2 


22.27 


19.33+0.06 


2.04+0.31 


0.77+0.09 


J111908.6+211917 


V 


15.79 


16.85 


0.71 


24.2 


6.92 


14.15+0.00 


-0.08+0.00 


0.35+0.00 


J112520.7+420425 


I 


19.36 


20.67 


0.28 


24.2 


21.47 


19.54+0.06 


1.76+0.45 


1.21+0.14 


J114526.3+193301 


V 


24.40 


22.67 


0.84 


24.8 


12.83 


22.79+0.36 


0.30+0.66 


0.88+0.65 


J114910. 5-002313 


V 


22.82 


23.44 


0.75 


23.3 


20.46 


22.62+0.95 


3.97+1.13 


0.41+1.30 


J114928. 3+143844 


V 


22.83 


20.71 


1.15 


24.5 


25.49 








J115642.8+022451 


V 


18.99 


20.42 


1.00 






18.21+0.02 


1.09+0.03 


0.32+0.03 




I 


17.81 


18.02 


1.00 


24.3 


10.35 


18.21+0.02 


1.09+0.03 


0.32+0.03 


J120326. 4+443635 


V 






1.43 


25.0 


0.0 








J121026.6+392909 


V 


19.53 


16.29 


1.27 


25.0 


12.55 


18.98+0.02 


0.34+0.03 


0.45+0.03 


J121658.4+375439 


I 


20.56 


19.36 


0.41 


23.3 


20.20 








J121705.5-031137 


I 


23.89 


22.49 


1.51 


24.3 


11.34 








J121707.7-031127 


I 


15.67 


17.34 


0.97 


23.2 


2.44 








J121839.7+295325 


V 






1.53 


24.8 


0.0 








J122331. 0+155245 


I 


21.02 


19.03 


0.64 


24.9 


18.67 


23.91+2.84 


0.92+0.77 


-1.66+2.87 


J122624.4+173228 


V 


18.82 


20.38 


1.15 


24.5 


13.56 








J125029. 2+302527 


I 


17.83 


19.93 


0.44 


23.9 


31.38 


17.85+0.02 


1.34+0.04 


0.53+0.03 


J125635.3+215632 


I 






0.52 


24.2 


0.0 








J125650.0+220630 


V 


22.44 


20.71 


0.83 


24.7 


5.71 








J131223. 6+424517 


I 


23.19 


22.90 


0.47 


25.9 


8.08 








J131617.8+420239 


V 


16.66 


19.42 


0.00 


24.3 


14.92 


16.12+0.01 


0.78+0.01 


0.45+0.01 


J134219.9-005816 


I 


24.08 


23.53 


1.23 


24.0 


17.42 








J140019. 9+044421 


I 


19.51 


20.32 


0.05 


24.6 


9.73 


19.51+0.10 


1.47+0.43 


1.08+0.22 


J143530.0+484534 


V 


22.07 


23.13 


0.94 


24.3 


17.58 








J150109.7+225111 


I 






0.82 


23.6 


0.0 








J151057.4+312722 


V 






0.49 


23.7 


0.0 








J153955.0+342013 


V 


19.49 


21.28 


1.11 


23.7 


26.62 


18.21+0.04 


1.86+0.12 


0.58+0.06 


J155938.7+473309 


V 


25.37 


24.08 


0.49 


25.0 


6.78 








J163141.4+373603 


V 


19.87 


20.87 


0.41 


24.7 


6.43 


18.57+0.05 


1.62+0.17 


0.81+0.08 


J163233.8+190550 


I 


22.36 


21.26 


0.26 


23.4 


17.40 








J164911. 4+305226 


V 






0.44 


24.3 


0.0 








J171901. 1+480458 


V 


21.01 


21.27 


0.34 














I 


19.57 


19.57 


0.34 


23.9 


11.65 








J172025.4+480321 


V 


21.66 


22.43 


1.48 


23.1 


75.08 








J172232. 9+501232 


V 


20.11 


20.65 


0.45 














I 


18.56 


19.44 


0.45 


24.4 


9.80 









a The name of the source as per IAU recommendations, as used 
by FIRST. FIRST Jhhmmss.s+ddmmss in which the coordinates 
are equinox J2000.0 and are truncated (Becker et al. 1995). 

h V and / refer to the HST-WFPC2 bandpasses F606W or F814W, 
respectively. 

c Surface brightness in mag/arcsec 2 within the effective radius 
(Re)- ' 

d Integrated flux at 1.4 GHz in mJy provided by the FIRST cat- 
alog. 

e The 50% completeness limit as determined by the Snigula et al. 
(2002) method (see § 2.4). 

f The likelihood ratio as discussed § 2.4. 

g Obtained from SDSS-dr5 (Adelman-McCarthy et al. 2007). 
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Fig. 2. — WFPC2 images of our 60 millijansky radio sources. Objects not located in the middle of the radio contours are indicated by 
an arrow. In the upper left of each stamp, we list the name from the FIRST catalog, the HST-WFPC2 optical magnitude in the given 
bandpass, the visually-defined morphology, and shifts in RA and Dec (Aa and A5) defined in § 2.3. Each stamp is 30" X 30 /x . Higher 
resolution images are available by request. 
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classification is listed for each object in Table 5. In order 
to also provide some quantitative measure of the galaxy 
morphologies, we fit the surface-brightness profiles with 
the Sersic or r 1 ^ 71 of model: 



I(r) = I e exp -b n 



l/n 



1 



(6) 



where I e is the intensity at the effective radius r e and b n 
is a constant defined in terms of n which describes the 
shape of the profile (eg. Graham & Driver 2005). The 
Sersic index (n) can morphologically classify galaxies by 
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type, where elliptical galaxies have n ~ 4 and later-type 
disk galaxies have n ~ 1. Therefore, we have a quanti- 
tative, morphological measurements to complement our 
visual classifications for the 25/43 classifiable, optical 
counterparts. We give these results in Table 5. 

For the 25/43 galaxies which are sufficiently bright and 
extended to be visually classified, there are 15 ellipticals, 
1 SO, 2 late-type spirals, 3 optical point sources (these are 
generally quasars, see § 4), and 4 mergers. Therefore, 
at the millijansky flux levels, 19/25 galaxies are early- 
type or quasars, while only 2/25 are later-type galaxies. 
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However, the surface brightness analysis permits a "con- 
tinuum" of galaxy types based on the best-fit value for 
the Sersic index. Therefore, to produce a comparable 
classification scheme (such as early- versus late- type), 
we collapse this continuum using a threshold value for 
the Sersic index of n > 2 are early-types and n < 2 are 
late- types (eg. Driver et al. 2006). Of the 34/43 galax- 
ies for which we could reliably measure the Sersic index, 
20 are early-types while 14 are later-types, and seven of 
these only have n<2 at the la-level (such as J002219.1- 
013030). Based on these Sersic indices, we find ~60% of 
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millijansky radio galaxies are morphologically early- type 
systems. Furthermore, these classifications based on the 
Sersic index generally agree with our visually assigned 
morphologies. Based on these morphological analyses, 
the dominance of red galaxies in millijansky radio galaxy 
samples (Kron, Koo, & Windhorst 1985; Hammer et al. 
1995a) is now confirmed by robust morphological obser- 
vations with high-resolution imaging from HST. 
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TABLE 5 

Morphological and Structural Results 



Radio Source a 


Re h 


n c 


2d 


ANN e 


Visual 1 


WFPC2S 




(arcsec) 






Type 


Type 


Target 


J002219.1-013030 


0.23 


1.56+0.82 


0.30 


2.5 


X 


ANY 


J004322. 3-001343 


0.11 






-4.1 


ps 


ANY 


J012616. 6-012126 


1.33 


4.95+0.74 


0.11 


5.0 


SBb-SBc 


ANY 


J030237.6+000818 










UNID 


FIELD-030239+00065 


J030249.5+000615 


0.51 


4.39+0.41 


0.01 


-2.3 


E 


FIELD-030251+00071 


J082820. 6+344321 


0.84 






5.5 


X 


6C0825+34 


J082828.4+344131 










UNID 


6C0825+34 


J084715.5+443752 


0.21 


1.00+0.00 


1.84 


8.5 


xE 


GAL-084720+443739 


J084849.5+445550 


0.90 


3.98+0.03 


0.19 


0.6 


E (cluster) 


LYNX-E 


J091205.2+350506 


0.22 


1.24+0.75 


0.05 


7.6 


xE 


ANY 


J091251.0+525928 


0.34 


1.80+0.70 


0.18 


5.2 


xE 


SBS0909+523 


J094926.6+295941 


0.91 


3.46+0.46 


0.01 


-3.1 


E 


ANY 


J094930.7+295938 


0.37 


2.00+0.54 


0.02 


1.1 


SO 


ANY 


J100354.5+285911 


0.86 






7.2 


merger 


ANY 


J102437.2+470312 


0.21 


3.77+0.02 


1.19 


-4.1 


xE 


PAR 


J102744.6+282921 


0.94 






6.8 


X 


HIGH 


J103452.3+394704 


0.44 


2.66+0.67 


0.20 


2.1 


xE 


PAR 


J104630.8-001213 


1.57 






3.0 


merger 


10HR-B 


J104757.0+123835 


0.43 


2.58+0.41 


0.01 


-2.5 


E 


ANY 


J111908.6+211917 


0.65 


2.63+0.23 


0.11 


-0.9 


ps 


PG1116+215 


J112520.7+420425 


0.73 


1.72+0.44 


0.01 


3.0 


E/S0 


HI-LAT 


J114526.3+193301 


0.18 


0.73+0.24 


1.08 


3.0 


xC 


ANY 


J114910. 5-002313 


0.53 


3.66+0.04 


0.50 


8.3 


xC 


ANY 


J114928. 3+143844 


0.15 


2.17+0.30 


0.25 


0.9 


E 


ANY 


J115642.8+022451 


0.44 


1.76+0.60 


0.04 


1.9 


E 


PAR 


J120326.4+443635 










UNID 


ANY 


J121026.6+392909 


0.09 


4.28+0.02 


0.73 


-2.8 


ps 


ANY 


J121658.4+375439 


0.23 


5.45+0.33 


0.06 


-2.6 


xE 


MS1214.3+3811 


J121705.5-031137 


0.21 


0.50+0.17 


0.00 


6.4 


X 


2MASSW-J1217-03 


J121707.7-031127 


0.86 


2.62+0.39 


0.03 


-3.5 


E 


2MASSW-J1217-03 


J121839.7+295325 










UNID 


ANY 


J122331. 0+155245 


0.16 






-1.7 


xC 


PAR 


J122624.4+173228 


0.82 


3.08+0.39 


0.04 


-2.2 


E (companion) 


ANY 


J125029. 2+302527 


1.05 


1.96+0.48 


0.09 


0.4 


E (companions) 


HI-LAT 


J125635.3+215632 










UNID 


ANYWHERE 


J125650.0+220630 


0.18 


3.81+0.07 


1.88 


-1.2 


xC 


ANY 


J131223. 6+424517 


0.35 


1.54+0.70 


0.10 


3.1 


E 


SSA13E 


J131617.8+420239 


1.42 


1.48+0.34 


0.54 


0.2 


SBc 


ANY 


J134219.9-005816 


0.31 






8.8 


X 


HI-LAT 


J140019. 9+044421 


0.58 


5.43+0.33 


0.06 


-0.5 


E (companions) 


BIG2 


J143530.0+484534 


0.65 






6.4 


merger 


NGC5689 


J150109.7+225111 










UNID 


TVLM513-46546 


J151057.4+312722 










UNID 


WAS96 


J153955.0+342013 


0.91 


2.68+0.54 


0.01 


3.6 


E 


ANY 


J155938.7+473309 


0.22 


0.47+0.19 


0.86 


8.6 


xC 


ANY 


J163141.4+373603 


0.63 


2.81+0.41 


0.04 


-3.3 


E 


PAR 


J163233.8+190550 


0.24 


1.88+1.72 


0.03 


3.4 


xC 


2MASSW1632+1904 


J164911.4+305226 










UNID 


HIGH 


J171901. 1+480458 


0.40 


3.10+0.51 


0.01 


-0.8 


E 


HIGH 


J172025.4+480321 


0.57 


0.68+0.28 


0.26 


6.1 


xE 


ANY 


J172232. 9+501232 


0.60 






-2.3 


merger 


PAR 



a The name of the source as per IAU recommendations, as used 
by FIRST. FIRST Jhhmmss.s+ddmmss in which the coordinates 
are equinox J2000.0 and are truncated (Becker et al. 1995). 

b Half light radius in arcsec as measured by LMORPHO (Odewahn 
et al. 1996). 

c Sersic index (n) defined by the equation (6). 

d Reduced xt between the observed profile and equation (6). 

e Artificial Neural Network type classification as measured by 
LMORPHO (Odewahn et al. 1996) 

f Morphology of optical counterpart, 'x' used to classify objects 
too small to assign a reliable morphology. 'xC' indicates the object 
is compact. 'xE' indicates the object is extended. '*' indicates a 
unique note for that source. 

g This is the header keyword for the field target from the WFPC2 
images, which helped us assess the "randomness" of each WFPC2 
field. 
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3.2. Optical Magnitude Distribution 

We measure optical fluxes from the HST-WFPC2 im- 
ages, using the software package LMORPHO (Odewahn 
et al. 1996), which we list in Table 4. In Figure 4, we show 
the distributions of the V- and /-band apparent magni- 
tudes. While the statistics may be relatively low, these 
distributions clearly peak at V^22 mag and J~26 mag, 
which is in each plot several magnitudes brighter than the 
typical field depth. These peaks are similar to the what 
was observed in the LBDS (Windhorst et al. 1984b), and 
its Hercules subfield (Waddington et al. 2000), which are 



only complete to V~22.7 mag and 7c^21 mag, respec- 
tively. Since our compendium of HST-WFPC2 fields are 
generally ~ 2 mag deeper than the optical imaging in the 
LBDS, these peaks in the magnitude distributions are 
more secure than reported in those earlier works, which 
were based on ground-based photographic or CCD imag- 
ing and poorer seeing. Furthermore, in the microjan- 
sky flux range, Hammer et al. (1995a) find similar peaks 
in the apparent magnitude distributions in the Canada- 
France Redshift Survey (CFRS). 

Since for a given object only HST/WFPC2 V- or I- 



12 



Russell et al. 



J1. 21 839+295325 
v Hiag 




8C5921 
AS* 0.4620 



■■.::=::: : -E:. \ ■ 




J1 22331 +155245 

L— &L0 mag ■*■■■ 
'.■ lAicla^ltiablei ■'. 
■^a=0;5921^ 
A8-- 0.4620 : 



Si# 



IHi; 



o: 




J 122624+173228 
: :; V,= 1 8.8 nriag : 
elliptical 

'■■};■;■■';':=' ;^a=;i. ■ 



^3 




J125p:29+302527 
*7.8gpg 



£ 



Aa=-U793 
A5=0,t394 



J1 25635+21 5632 
I = ... mag 
Unid 

Ao=- 1.0356 
Ap=0,0707 







band imaging is available, we require an algorithm to 
convert the observed flux in one bandpass to what is 
expected in the other. Therefore, we retrieve the </-, 
r'-, and z'-fluxes for 22/43 galaxies available from the 
Sloan Digital Sky Survey, Data Release 5 (SDSS-DR5; 
Adelman-McCarthy et al. 2007). The g'- and r'-band 
fluxes are converted to an effective F-band for a typi- 
cal galaxy spectral energy distribution (SED) following 
Jester et al. (2005): 

Kff =g'~ 0.59 x (g f - r') - 0.01 (mag). (7) 

In Figure 5, we show the SDSS i'-band flux as a func- 
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tion of the effective F-band flux, as synthesized from the 
SDSS g'- and r'-band observations and the best-fit, linear 
model to fit of i' = (1.85 ± 0.02) + (0.86 ± 0.00) x F eff . To 
verify these conversions from g' and r' to V e R and from 
i' to /, we show AV = (V - V eS ) and Al = (I - i') as 
a function the respective HST observations in Figure 6. 
For galaxies which only have V-band WFPC2 data, the 
difference between the HST-WFPC2 and the effective V 
is (AV) = 0.03 ±0.18 mag, while the corresponding differ- 
ence in the /-band is (Al) =0.03 ± 0.25 mag. This indi- 
cates that the correlation in Figure 5 can be used to sys- 
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tematically convert the WFPC2 F-band fluxes to corre- 
sponding /-band fluxes with an uncertainty of < 0.3 mag. 
This uncertainty is not large enough to affect any of our 
conclusions in a significant way. 

With these converted fluxes, we show the /-band mag- 
nitude distribution for the 43 identified optical galaxies 
in Figure 7. To determine the location of the peak, we fit 
a simple Gaussian to magnitude distribution. The peak 
at 7:^20.7 ± 0.5 mag is significantly brighter (~3 mag) 
than the completeness limit of our typical WFPC2 field, 
which suggests that this peak is not due to incomplete- 



ness. This peak may reflect the redshift distribution, 
and possibly strong cosmological evolution of the radio 
source population. Waddington et al. (2000) point out 
that if the radio galaxies had the same space density at 
all redshift s, then we expect the number of galaxies per 
magnitude interval would increase toward the observa- 
tional limit, as is roughly the case for the optical field 
galaxy counts. With the exceptional resolution and in- 
creased depth of the HST-WFPC2 dataset, we extend 
the typical limit of AB ~ 22 mag to AB ~ 24 mag for 
the millijansky population over a significant field-of-view. 
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This magnitude distribution is quite different than what 
is expected from a general population of field galaxies 
(Tyson 1988). The radio-selected galaxy counts likely 
peak at 7:^ 20.7 ±0.5 mag since the majority of millijan- 
sky sources have L*-type luminosities at z^O.S (Condon 
1989). 

3.3. Radio-to- Optical Spectral Index 

Kron, Koo, & Windhorst (1985) show that the radio- 
to-optical spectral index (a ro ) can be used to distinguish 
between red and blue galaxies. Therefore, we perform 



a similar analysis with morphological types determined 
from the high-resolution imaging. In Figure 8, we show 
the /-band flux as a function of the radio flux for the 
43 optically identified galaxies. Objects are separated 
into three classes by their optical morphology (discussed 
in § 3.1) of early, late, and unclassifiable as filled red 
circles, green triangles, and blue squares, respectively. 
We show lines of constant radio-to-optical spectral index 
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Fig. 3. — Distribution of likelihood ratios, LR as defined by 
equation (1). We compute the likelihood that a radio source has a 
given optical counterpart following de Ruiter et al. (1977), which 
is summarized in § 2.4. We show all radio sources which were posi- 
tively identified with an optical counterpart as a solid distribution. 
The unidentified galaxies are represented as a hatched distribution. 
As argued in the text, every positively identified object has LR>2, 
while the unidentified galaxies have LR<^2. Note that two objects 
with LR > 40 are not represented. 

defined by: 




ol to — 



log (Si. 4 /Sj) 



log (20.4 cm/8012 A) 
0.185 x log (S! A / Si), 



(8) 
(9) 



where log (Sj/Uy) = -0.4 x (7-8.72) mag (Oke & Gunn 
1983). The solid lines denote logarithmic radio-to-optical 
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Fig. 4. — Optical WFPC2 magnitude distributions, separated by 
bandpass (top: V-band, bottom: /-band). The arrows represent 
the depth for fields where no optical counterpart could be identified. 
The vertical dashed line shows the flux- averaged, 50% complete- 
ness limit for each bandpass. These magnitude distributions peak 
>2 mag brighter than the average completeness limit for the aver- 
age WFPC2 field. While the statistics may be relatively low, this 
suggests that the turn-over in the optical magnitude distributions 
from earlier studies was not due to their shallow imaging. 

flux ratios from 1 to 10 5 . In general, the elliptical galax- 
ies (filled red circles) have / < 20 mag and shallower 
spectral indices of a ro < 0.4. Conversely, the unclassifi- 
able galaxies (filled blue squares) have a ro > 0.4, which 
is primarily due to their brightnesses (7>21 mag). Fur- 
thermore, since these galaxies span the entire range of 
radio fluxes, they are likely representative of the gen- 
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v eff (mag) 
Fig. 5. — Colors of FIRST radio galaxies in our sample from 
SDSS-DR5. The effective V-band is constructed following Jester 
et al. (2005) as a combination of the SDSS g' - and r'-bandpasses, 
which we give in the lower-right corner of the upper panel. The 
solid line represents the best-fit linear model to the data, which is 
given in the upper-left corner of the upper panel. The lower panel 
shows the residuals, as the observations minus the model. The 
average color for these radio galaxies from the SDSS database is 
((i 7 - Veff)) -1.8 mag. 
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Fig. 6. Comparison of the fluxes from SDSS-DR5 and HST- 
WFPC2 observations. These panels show the ratio between the 
HST and the SDSS fluxes as a function of the HST flux. The 
average difference and standard deviation are listed in the lower- 
right of each panel. In the left panel, we show the (F606W — V e n) 
residual for the objects which were only observed in the V-band 
and had SDSS data. Since the average difference is only —0.08 ± 
0.79 mag, the conversion from g' and r' to V^ff, and hence the 
correlation shown in Figure 5 is relatively robust. The right panel 
shows that the once the V-band measurements are converted to 
the i'-band, there is relatively little additional error introduced by 
assuming that i 1 ~I for these radio galaxies. 

eral radio sample. Their lower optical fluxes could be 
caused by an increased dust obscuration associated with 
a surrounding starburst and/or enhanced/beamed radio 
emission. Deep infrared observations would be capable 
distinguishing these two scenarios, and likely detect the 
8/43 optically unidentified radio galaxies. 

In the GOODS-S field (Giavalisco et al. 2004), Afonso 
et al. (2006) publish the optical and radio brightnesses for 
their sample of 43 optically-identified micro jansky radio 
galaxies with morphological classifications. In Figure 8, 
we show their micro jansky radio galaxies as open sym- 
bols, where late- type and star- forming galaxies generally 
constitute a much larger fraction of radio galaxies than 
at the millijansky level. Moreover, these late- type galax- 
ies generally have a ro < 0.2, which is a region of (a ro — I) 
parameter space where few early-type galaxies are found 
in our millijansky sample. For galaxies sufficiently re- 
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Fig. 7. — /-band magnitude distribution for all identified radio 
galaxies. The fluxes of sources where only V-band observations are 
available are converted assuming the SDSS-DR5 colors described 
in § 3.2, Figure 5, and Figure 6. The vertical dashed line shows 
the flux-averaged, 50% completeness depth of the the fields, while 
the arrows represent the depths of the eight fields where no optical 
counterpart was identified to their respective field limits. This peak 
at /~ 20. 7 ±0.5 mag is nearly 3 mag brighter than the depth of the 
typical WFPC2 field. The source counts of a random population of 
field galaxies should increase as a power-law to the completeness 
limit of the field, contrary to what is seen for the radio selected 
galaxies. 



SO, ellipticals, QSOs • 
spirals * 
unclassiriable ■ 
mergers * 




1 

log S 14 (mJy) 

Fig. 8. — Integrated radio flux versus total /-band magnitude for 
the complete FIRST+WFPC2 sample. The V-band flux has been 
converted to an effective /-band as described in § 3.2. Objects are 
separated into early, late, merger, and unclassifiable types, while 
the filled and open symbols are from this work and Afonso et al. 
(2006), respectively. The lines are contours of constant radio-to- 
optical flux ratios as defined by equation (8). In the millijansky 
flux range, the majority of the galaxies have early-type morpholo- 
gies and shallow spectral slopes of a ro ^0.4. Additionally, there is 
only one late-type spiral and 4 mergers which span a large range 
of spectral indices (0 < a ro < 0.5). However, lower radio fluxes 
(Si. 4 <C 1 mJy), there are very few early- type systems where late- 
type systems dominate the counts with the a ro ^0-4. This dra- 
matic change in optical morphologies at Si .4 ~ 1 mJy may be 
related to a change in the source of the radio emission, such as 
primarily starburst galaxies (Si .4 <1 mJy) versus primarily AGN 
(Si a ^1 mJy). Additionally, the line of a ro — 0.18 marks a fiducial 
point, below which nearly all galaxies are late- type. 

solved and given a reliable optical classification, spirals 
and mergers dominate the radio source counts at the mi- 
crojansky level (25/43, 53%) while early-type galaxies 
are the dominant type at millijansky fluxes (19/25, 76%). 

3.4. Optically Unidentified Radio Sources (Unids) 

For the sample of 51 FIRST sources, eight were found 
to have no optical counterpart identifiable to AB ~ 
23.6 - 25.0 mag (see Table 4 for 50% depths of each 
field from the Snigula method described in § 2.4). Like 
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Fig. 9. — The Hubble diagram for faint radio galaxies. Data is 
compiled from numerous sources (Windhorst et al. 1984b; Kron, 
Koo, & Windhorst 1985; Windhorst et al. 1985; Fomalont et al. 
1997; Richards et al. 1998; Waddington et al. 2001; Roche et al. 
2002; Chapman et al. 2003). The circles, triangles, and squares 
represent ellipticals, starbursts, and quasars, respectively. The 
open symbols represent point-like or small radio sources, while 
filled symbols show (classical) double or extended radio sources; 
the size of the plotted point is logarithmically proportional to the 
radio flux from 10 — 10 5 microjansky (Windhorst 2003). The lines 
are from stellar population synthesis models of (Bruzual &; Char- 
lot 2003). The dot-dashed line is a constant star-formation model 
and solid, dashed, and dotted lines are an exponential star-forming 
history with ages of 13, 11, and 9 Gyr, respectively. The e-folding 
timescale in the exponential models is 1.3 Gyr. All models assume 
a stellar mass of 10 11 Mq, a Salpeter IMF (Salpeter 1955), and 
solar metallicity. The obvious trend in this diagram suggests that 
early- type radio galaxies are roughly standard candles, and that 
the unidentified millijansky objects with / > 22 mag are likely at 
z>l. NOTE: The /-band magnitudes in this figure are in the Vega 
system, when necessary we assume (Jab — ^Vega) — 0-48 mag. 

the unclassifiable galaxies, the unidentified objects have 
radio fluxes which span a comparable range to the entire 
radio sample (see Table 3) and thus may be represen- 
tative of the general population. However, unlike the 
unclassifiables, their lack of any optical flux may indi- 
cate that they are are considerably higher redshift. In 
Figure 9, we show the Hubble diagram (/-band flux as 
a function of redshift) for all faint radio galaxies, where 
the size of the plotted symbol is logarithmically propor- 
tional to the radio flux (see Windhorst 2003, for more 
details). From Figure 9, millijansky radio sources with 
an optical counterpart with a brightness of / > 22 mag 
are likely at z > 1 (eg. Schmidt, Connolly, & Hopkins 
2006), indicating that these galaxies may be standard 
candles. However, the significant spread in redshift for 
galaxies with i>22 mag suggests than an additional ef- 
fect may be present, which could be an increased optical 
extinction, different stellar populations at high redshift, 
or significant luminosity evolution. Without a robust 
optical detection of these galaxies, we can only speculate 
about their counterparts. Based on the Hubble diagram, 
these optically unidentified galaxies could be similar to 
the distant red galaxies (DRGs; Forster Schreiber et al. 
2004; van Dokkum et al. 2006; Georgakakis et al. 2006; 
Schmidt, Connolly, & Hopkins 2006) at z~2. The DRGs 
are generally old, dusty, and massive galaxies with very 
red infrared colors {I — J)>2 mag, ( J — H) ~ 2 mag, and 
(H — K) ~ 1 mag. Therefore, near-infrared observations 
are ideal for identifying the optical/infrared counterpart. 



4. COMMENTS ON INDIVIDUAL RADIO SOURCES 

We give brief comments on each of the 51 radio sources 
in the sample and their most likely optical identification 
(where applicable). We note the overlap with previous 
works where available found from the NASA/IPAC Ex- 
tragalactic Database 4 . 

1. J002219. 1-013030: Unclassifiable. The optical 
counterpart is entirely visible despite being near 
the edge of the WFPC2 field. The object is also 
within 18" of a NVSS radio source (Condon et al. 
1998). 

2. J004322. 3-001343: Point-source. This object 
was identified as a quasar in the SDSS (Richards 
et al. 2004). 

3. J012616. 6-012126: Late-type, barred spiral 
(type SBb-SBc). This object is a bright (V ~ 
19.8 mag) and well-resolved mid-late type spiral 
galaxy. 

4. J030237.6+000818: Unidentified with I > 
24.5 mag. There are four CFRS objects within 30" 
(Hammer et al. 1995b). 

5. J030249. 5+000615: Elliptical. There is no no- 
ticeable optical elongation to indicate any align- 
ment of the radio source within the galaxy. Ham- 
mer et al. (1995b) classified this source as an E0. 

6. J082820. 6+344321: Unclassifiable. The optical 
candidate is very faint (/c^25.6 mag). 



7. J082828.4+344131: 
24.7 mag. 



Unidentified with / > 



8. J084715. 5+443752: Unclassifiable. The radio 
source is clearly extended and was found by Wind- 
horst et al. (1985). 

9. J084849.5+445550: Elliptical. There is no evi- 
dence for any optical elongation with the resolved 
radio source. The object is located in a group or 
poor cluster, and the shape of radio contours were 
used to determine the identification. This object 
was also identified by Oort (1987). 

10. J091205. 2+350506: Unclassifiable. The optical 
source is extended and very faint (V~23.8 mag), 
and is near a bright star. 

11. J091251. 0+525928: Unclassifiable. The optical 
counterpart is small, and possibly extended. The 
radio counters are clearly elongated. 

12. J094926. 6+295941: Elliptical. The optical and 
radio emission is not aligned. 

13. J094930.7+295938: SO. The radio source is 
clearly double-lobed and the optical identification 
is slightly off the radio axis, the radio contours ap- 
pear to be somewhat "bent" toward the identified 
object. 

4 NED is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aero- 
nautics and Space Administration. 
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14. J100354.5+285911: Merger. The radio source 
is likely double-lobed with the merging system lo- 
cated at the center of the radio axis. 

15. J102437.2+470312: Unclassifiable. The optical 
object is marginally extended. 

16. J102744.6+282921: Unclassifiable. This object 
is very faint (J c^ 22.4 mag) and on the outskirts 
of a bright spiral galaxy, which is ~ 7" to the 
west. Therefore, the spiral galaxy center is likely 
not the optical identification. The radio source may 
be caused by a supernova remnant (eg. Windhorst 
et al. 1984b; Kron, Koo, & Windhorst 1985) or a 
background galaxy, which may be somewhat ex- 
tincted by the foreground spiral galaxy. 

17. J103452. 3+394704: Unclassifiable. The optical 
source is extended, likely a galaxy in a group. 

18. J104630. 8-001213: Merger. The merger pair is 
an elliptical and early-type disk galaxy, however 
neither show obvious signs of interaction. The op- 
tical candidate was in the sample of Glazebrook et 
al. (1995). 

19. J104757.0+123835: Elliptical. The radio source 
is fairly round and was identified by Condon et al. 

(1998). 

20. J111908. 6+211917: Point source. This object 
was identified as a QSO by Barkhouse & Hall 
(2001). 

21. J112520. 7+420425: Elliptical/SO. The radio 
source is round and was identified by Condon et 
al. (1998). 

22. J114526. 3+193301: Unclassifiable. 

23. J114910. 5-002313: Unclassifiable. 

24. J114928. 3+143844: Elliptical. This galaxy was 
identified by Condon et al. (1998). 

25. J115642. 8+022451: Elliptical. The optical emis- 
sion is roughly perpendicular to the radio contours. 
This object was identified by Condon et al. (1998). 

26. J120326.4+443635: Unidentified with V > 
25.0 mag. 

27. J121026. 6+392909: Point source. It is located 
in group and was classified as a QSO (Hewitt & 
Burbidge 1993). 

28. J121658.4+375439: Unclassifiable. The optical 
source is marginally extended. The radio source 
may be double-lobed, but of very unequal lobe flux. 

29. J121705. 5-031137: Unclassifiable. The radio 
source is resolved. 

30. J121707.7-031127: Elliptical. The optical emis- 
sion is round, while the radio counters are elon- 
gated. This object was classified as an AGN in the 
SDSS (Richards et al. 2004). 



31. J121839.7+295325: Unidentified with V > 
24.8 mag. The morphology of the nearby opti- 
cal source suggests that it could be gravitationally 
lensed by the unidentified radio galaxy (Ryan et al. 
2008b). 

32. J122331. 0+155245: Unclassifiable. This radio 
source was identified by Condon et al. (1998). 

33. J122624.4+173228: Elliptical. This galaxy may 
have a companion to the east-north-east, and is in 
a small group. This radio galaxy was identified by 
Condon et al. (1998). 

34. J125029. 2+302527: Elliptical. This galaxy is 
like in a cluster or group. The radio source has 
round inner contours, but appears to be resolved 
on larger scales. 

35. J125635. 3+215632: Unidentified with / > 
24.2 mag. There are three relatively bright objects 
within ~ 4" of the radio source. The radio source 
was identified by Condon et al. (1998). 

36. J125650. 0+220630: Unclassifiable. The radio 
source was identified by Condon et al. (1998). 

37. J131223.6+424517: Elliptical. 

38. J131617.8+420239: Late-type barred spiral 
(SBc). This object is >10" in diameter and is also 
an infrared source, identified in the Two Micron All 
Sky Survey (2MASS; Skrutskie et al. 2006). 

39. J134219. 9-005816: Unclassifiable. 

40. J140019. 9+044421: Elliptical. The object is lo- 
cated in the center of a small group. The radio 
contours are slightly elongated, while the optical 
emission is fairly round. 

41. J143530. 0+484534: Merger. The optical candi- 
date is on innermost radio contour. 

42. J150109. 7+225111: Unidentified with / > 
23.6 mag. There is a disk galaxy c^3" from the ra- 
dio position that could be the optical counterpart. 
Since the astrometric correction for this field was 
small, the nearby disk galaxy is likely not related, 
owing to the large separation and the relatively low 
background object density. 



43. J151057.4+312722: 
23.7 mag. 



Unidentified with V > 



44. J153955. 0+342013: Elliptical. The optical emis- 
sion is fairly elongated, while the radio contours are 
round. This system may have several companion 
systems to the south and west. 

45. J155938. 7+473309: Unclassifiable. This object 
may be located in a dense group of faint, small ob- 
jects. This radio source was identified by Condon 
et al. (1998). 

46. J163141.4+373603: Elliptical. The optical and 
radio emission appear to be aligned. The radio 
emission is part of a larger, low SB complex. 
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47. J163233. 8+190550: Unclassifiable. The optical 
counterpart is the northern of the two small ob- 
jects. This radio source was identified by Condon 
et al. (1998). 

48. J164911.4+305226: Unidentified with V > 
24.3 mag. There is a late type galaxy approxi- 
mately 4" from the radio position, but given that 
the astrometric corrections are less than 2", and the 
LR ~ 0.0, this is considered an unidentified source. 

49. J171901. 1+480458: Elliptical. This radio source 
found in Condon et al. (1998). 

50. J172025.4+480331: Unclassifiable. This object 
may be marginally extended. 

51. J172232. 9+501232: Merger. The primary 
galaxy is clearly a late-type, possibly barred spi- 
ral which shows clear signs of interaction. The sec- 
ondary galaxy is considerably smaller and may also 
be a late-type disk system. 

5. SUMMARY AND DISCUSSION 

We have determined the optical morphology and struc- 
ture from HST-WFPC2 imaging of a randomly-selected 
sample of 51 millijansky radio sources from the FIRST 
survey. Earlier studies of millijansky radio galaxies have 
often distinguished early- versus late-types based on ob- 
served optical colors or spectra. However, we have ex- 
tended such works by providing optical classifications 
based on HST imaging, and confirm that elliptical galax- 
ies constitute the majority of the millijansky source pop- 
ulation. 

We find that the optical flux distribution of these 
galaxies peaks at 7^20.7 ± 0.5 mag, which is consistent 



with evolutionary model of Condon (1989): faint radio 
galaxies generally have L*-type optical luminosities with 
a median redshift of z~ 0.8 and an effective redshift de- 
cline for z > 1.5. At low redshift (z < 1), most radio 
sources brighter than S1.4 > 1 mJy are generally found 
in massive ellipticals, and there is a clear deficit of these 
massive galaxies at z>3 (eg. Driver et al. 1998). De Lu- 
cia et al. (2006) argue that ~50% of the stars which will 
likely end up in elliptical galaxies in the local Universe 
are formed at z c^ 3, and are not finally assembled into the 
galaxy until zc^O.S. This down-sizing picture requires an 
increased major merger rate from z ~ to z ~ 3, which 
has been observed in the Hubble Ultra Deep Field (eg. 
Ryan et al. 2008a). Finally, this effective redshift cut-off 
may arise, since ~ 1 Gyr may be on average required to 
trigger an active nucleus after the major merger based 
on hydrodynamical models (Springel et al. 2005). This 
must be further studied in the critical redshift range of 
1 < z < 3, which may prove useful in tracing the high- 
redshift evolution of the supermassive black hole-bulge 
mass relation (eg. Ferrarese & Ford 2005). At these red- 
shifts, the optical counterparts will be extremely faint 
(AB > 25 mag), and may have a high dust content or their 
luminosity function may evolves strongly or episodically 
with redshift. Consequently, dedicated observations with 
the Wide Field Camera 3 (WFC3) or the James Webb 
Space Telescope will be needed for the continued study 
of the population of millijansky radio sources. 
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